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In this Letter we consider a liquid mixture confined between two thermally conducting walls 
subjected to a stationary temperature gradient. While in a one-component liquid non-equilibrium 
fluctuation forces appear inside the liquid layer only, non-equilibrium fluctuations in a mixture 
induce a Casimir-like force on the walls. The physical reason is that the temperature gradient 
induces large concentration fluctuations through the Soret effect. Unlike temperature fluctuations, 
non-equilibrium concentration fluctuations are also present near a perfectly thermally conducting 
wall. The magnitude of the fluctuation-induced Casimir force is proportional to the square of the 
Soret coefficient and is related to the concentration dependence of the heat and volume of mixing. 
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When large and long-range fluctuations are present, 
they will induce forces in confined fluids [T]. These 
are commonly referred to as Casimir-like forces in anal¬ 
ogy to forces induced by vacuum fluctuations between 
two conducting plates mm- A well-known example is 
the Casimir force induced by critical fluctuations in flu¬ 
ids [HZ]- Apart from critical systems, long-range corre¬ 
lations also exist in equilibrium systems with Goldstone 
modes ^ and in many non-equilibrium systems, where 
even longer-range correlations can exist [mu. 

In this Letter we consider a liquid mixture in a non¬ 
equilibrium steady state (NESS) between two parallel 
thermally conducting plates subjected to a uniform tem¬ 
perature gradient VT. In a liquid mixture a temper¬ 
ature gradient induces large concentration fluctuations 
through the Soret effect min]. These non-equilibrium 
concentration fluctuations vary with the 4*** power of 
the inverse of the wave number k of the fluctuations, 
just as the non-equilibrium temperature fluctuations in 
a one-component fluid |S1[T1]. However, there is a prin¬ 
cipal difference between the Casimir pressures induced 
by non-equilibrium concentration fluctuation and those 
induced by non-equilibrium temperature fluctuations. In 
thin fluid layers, fluctuations not only may induce a force 
on the walls, but also may introduce an effective po¬ 
tential inside the fluid layer causing a modification of 
the density or composition profile m- While in a one- 
component fluid non-equilibrium fluctuations only induce 
the latter phenomenon yielding a re-arrangement of the 
density profile dl’ the purpose of the present letter is to 
demonstrate that non-equilibrium concentration fluctua¬ 
tions induce an actual Casimir pressure on the walls. 

It is well known that in considering the dynamics 
of fluctuations around thermal equilibrium, nonlinear 
terms in the hydrodynamic equations serve to renormal¬ 
ize various terms in the linearized hydrodynamic equa¬ 
tions [17H2^ . Here we show that in a NESS the non¬ 
linear terms cause a most important renormalization of 


the non-equilibrium (NE) pressure or normal stresses in 
a binary fluid. To determine the non-equilibrium induced 
pressure in a liquid mixture, we need to consider the pres¬ 
sure p as a function of the fluctuating conserved quanti¬ 
ties, which are the fluctuating energy density e J- 5e, the 
fluctuating mass densities pi + Spi, and P 2 + Sp 2 of com¬ 
ponents 1 (solute) and 2 (solvent). As in the case of a 
one-component fluid, we can neglect the fast propagating 
sound modes and, hence, the linear fluctuation contribu¬ 
tion to the pressure m- Applying a Taylor expansion 
to the pressure then yields a contribution quadratic in 
terms of 6e, Spi, and Sp 2 ' 
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with a = (e, pi, p 2 )- In a liquid mixture there are two 
diffusion modes that are linear combinations of heat dif¬ 
fusion and mass diffusion naEg. An important param¬ 
eter for dealing with fluctuations in liquid mixtures is 
the Lewis number, which is the ratio of thermal diffu- 
sivity Dt and mutual mass diffusivity D-. Le = Dp/D. 
In liquid mixtures this Lewis number is commonly larger 
than unity. Hence, in dealing with fluctuations in liq¬ 
uid mixtures one often adopts a large-Lewis-number ap¬ 
proximation m- For large values of the Lewis number, 
these diffusion modes decouple into a pure temperature 
fluctuation mode with a decay time proportional to 
and a concentration fluctuation mode with a decay time 
proportional to D~^ [26]. Hence, to get the slowest fluc¬ 
tuation mode contribution for Le ^ 1, we not only may 
neglect linear pressure fluctuations, but also linear tem¬ 
perature fluctuations. 

For the concentration variable we adopt the mass frac¬ 
tion of the solute w = pxjp with p = pi + P 2 being 
the mass density of the mixture. At constant p and 
T, 5e, 6pi Sp 2 are related to the concentration fluctua- 
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tions 6 w by 5e = (de/dw)^ j, 6 w, Spi = {dpi/dw)^ 6 w, 
6 p 2 = {dp 2 /dw)pj,Sw. We then obtain from Eq. Q 
for the average NE contribution (r) at a position 
r = {x,y,z} to the equilibrium pressure p in terms of 
e,p,w: 


Pne(’') 


l/^\ (^\_f9e\ fd^p\ 

2Uejp,^ [\dwyp,T \9p)p,Sdw^)p,T 



where the superscript w indicates that (r) is a 
Casimir pressure induced by concentration fluctuations. 
We note that only the NE concentration fluctuations 
{{ 6 w{r)r) NE cause a renormalization of the pressure, 
since the equilibrium concentration fluctuations are al¬ 
ready incorporated in the unrenormalized pressure. Just 
as for the case of a one-component fluid [16], the NE 
pressure can be obtained from an explicit mode-coupling 
theory generalized to NESS, which justifies the approach 
adopted above. 

Relevant thermodynamic relations, associated with the 
hydrodynamic modes in a mixture, can be found in an 
article of Wood [3^. Noting that the thermodynamic 
field conjugate to the mass fraction w is the difference 
between the specific chemical potentials of the solute and 
the solvent, p = pi — p 2 , we can transform Eq. ([^ into 
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where a = —p~^(dp/dT)^^ is the thermal expansion 
coefficient, Cp^w the isobarlc specific heat capacity, 7 = 
cv,wlcp^w the ratio of the isochoric and isobaric heat ca¬ 
pacities, and X = idw/dp)p rp an osmotic susceptibility. 
This osmotic susceptibility can be related to the molar 
excess Gibbs energy m and, hence, its temperature and 
pressure derivatives to the excess molar enthalpy and 
the excess molar volume E®, so that 


Pne(^) 


P(7-1) 

2aT MfMi 
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where Mi and M 2 are the molar weights of the solute 
and solvent, respectively, M = MiXi -f (1 — xi)M 2 the 
molar weight of the mixture, and xi the mole fraction of 
the solute. 

The Intensity of the NE concentration fluctuations 
{{5w(r)f) NE can be obtained by solving appropriate fluc¬ 
tuating hydrodynamics equations uniiMiiis]. We con¬ 


sider a liquid mixture subjected to a stationary temper¬ 
ature gradient VT confined between two horizontal ther¬ 
mally conducting plates located at z = ±L /2 in the co¬ 
ordinate direction perpendicular to the plates. Such a 
temperature gradient Induces a stationary concentration 
gradient Vw = —w (1 — w) St'^T, where St is the Soret 
coefficient [niin]. A procedure for solving the fluctu¬ 
ating hydrodynamic equations to obtain the intensity of 
the NE concentration fluctuations has been developed 
by two of us, but with artificial boundary conditions for 
the fluctuations at the walls adopted for mathematical 
convenience [29] . It turns out that exactly the same pro¬ 
cedure can be used to obtain the solution for the inten¬ 
sity of the NE concentration fluctuations which satisfies 
physically realistic boundary conditions, namely, a rigid¬ 
boundary condition for the wall-normal velocity fluctu¬ 
ations, 6 vz = dSvz/dz = 0 1301, and the condition of 
no mass flux through the boundaries, d 6 w/dz = 0 , at 
2 ; = ±T/2, for Le S> 1. The solution (((5t(;(r))^)NE for 
the Intensity of the NE concentration fluctuations and, 
hence, for the NE Casimir pressure p^^ (z), only depends 
on the height z in the liquid layer. While the solution for 
arbitrary values of z is rather complicated, the important 
new result is that we have obtained a simple expression 
for the concentration fluctuations at the walls: 

{{Sw (z = ±^))\= (1 - StL i^Tf 
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where /cb is Boltzmann’s constant and u Is the kinematic 
viscosity. In Eq. (§ q — k\\L, where /cy is the magni¬ 
tude of the component of the wave vector k of the fluc¬ 
tuations parallel to the plates |33|. A derivation of the 
intensity of the NE concentration fluctuations, given by 
Eqs. and ([^, can be found in |3l|. Substitution of 
Eq. ([^ into Eq. @ yields for the NE Casimir pressure 
Pne — PPl‘P) exerted on the walls: 
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It has been verified experimentally that approximating 
the thermodynamic and transport properties In Eq. ([^ 
and, hence in Eq. Q , by their average values In the cen¬ 
ter of the liquid layer reproduces the Intensity of the non¬ 
equilibrium fluctuations to within 1 % at temperature dif¬ 
ferences up to AT = 40 K between two plates I 1 I 1 I 321 . 
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We note that for a given temperature gradient VT, the 
Casimir pressure exerted on the walls increases with the 
distance L between the plates, indicating that we are 
dealing with a giant, i.e., surprisingly large, Casimir ef¬ 
fect [53]. The physical reason is that the NE correla¬ 
tions diverge as which means that in real space the 
correlations scale with the system size L. While there 
exists an extensive literature on long-range correlation 
in NESS, we emphasize that only NE temperature and 
NE concentration fluctuations cause such a dramatic ef¬ 
fect liiin]. Physically, it may be more practical to study 
the NE Casimir pressure as a function of the distance L 
for a given temperature difference AT between the plates 
so that VT = AT/L. It then follows from Eq. Q that 
Pne be proportional to L~^ (AT/T)^. In principle 
Pne is also affected by gravity. However, it is expected 
that this effect will be minor except when the mixture is 
close to a hydrodynamic instability [T5] . 

It is interesting to compare the Casimir pressure in¬ 
duced by NE concentration fluctuations with a Casimir 
pressure induced by NE temperature fluctuations in a 
one-component fluid m- In analogy to Eq. Q, that 
result can be written as. 
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where h is the specific enthalpy and v the specific vol¬ 
ume. While Eqs. Q and Q look very similar, there is a 
fundamental difference between the two. In contrast to 
Eq. Q, Eq. implies 

Pne = ^"2 ) “ 


since at a thermally conducting wall 6T = 0. Hence, 
the Casimir pressure induced by NE temperature 
fluctuations only appears in the inside of the fluid layer, 
causing a re-arrangement of the density profile m, but 
it does not exert a pressure on the thermally conducting 
walls. Hence, we only quoted finite values for (pne)^ 
averaged over the height of the fluid layer in previous 
publications [lU [33] . We note that in the dilute-gas limit 
Pne vanishes for any value of z. 

To understand the implications of Eq. Q, we may en¬ 
vision a configuration where a (thin) plate with tempera¬ 
ture T 2 is located in a liquid mixture between two walls, 
both with a temperature Ti, as schematically illustrated 
in Fig. When > 0, the liquid mixture will ex¬ 
ert NE Casimir pressures on the two sides of the inner 
plate proportional to (AT)^/Ti and (AT)^/L 2 . When 
7 ^ L 2 , the plate will experience a net force causing it 
to move to the center of the liquid mixture layer. Hence, 
the force needed to move this plate off center, would be 


TABLE I: Estimated Casimir pressures 




L = 10"® m 

L = 10~^ m 

Pern 


-1 X 10"® Pa 

-1 X 10"“ Pa 

Pc* 


-6 X Pa 

-6 X 10"“ Pa 

Pne 

toluene-|-n-hexane“ 

+2 X 10"^ Pa 

+2 X 10“® Pa 

Pne 

l-methylnaphtalene-|- 

-1-9 Pa 

+0.9 X 10“^ Pa 


n-heptane'^ 



Pne 

aniline -I- methanoE 

-3 X 10"^ Pa 

-3 X 10"® Pa 


“Refs. 1^1161 
'’Ref. [TB] 

“Equimolar mixture, T = 298 K, AT = 25 K 


a measure of the Casimir force induced by the NE con¬ 
centration fluctuations. In practice it may be difficult 
to maintain plates at a close distance parallel to each 
other |3l|34l|35] and one may want to replace the plate 
by a particle. While a geometrical analysis of such a con¬ 
figuration becomes more complicated jSldS], the physical 
principle remains the same. 

In Table [l| we present some order-of-magnitude esti¬ 
mates for the pressure induced by NE concentration 
fluctuations and compare the values with those for the 
original Casimir pressure Pemj originating from electro¬ 
magnetic fluctuations in vacuum |5], and with the criti¬ 
cal Casimir pressure Pc mm- The Casimir pressures in¬ 
duced by NE concentration fluctuations are much larger, 
than either pem or pc. Table |l] also shows that pj^g can 
be either positive or negative, essentially depending on 
whether the concentration dependence of the heat of mix¬ 
ing is convex or concave. 



FIG. 1: Schematic illustration of Casimir pressures > 
0, induced by NE concentration fluctuations on a plate with 
temperature T 2 located in a liquid mixture between walls with 
temperature Ti. For < 0, the plate would be pulled to 
the closest wall. 
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Some authors have proposed NE Casimir forces in¬ 
duced by long-range fluctuations that originate from the 
spatial dependence of the noise correlations associated 
with the local fluctuation-dissipation theorem in the pres¬ 
ence of a gradient [3S1137]- However, for fluids it has 
been found that these NE correlations are insignificant 
compared to the NE correlations caused by a coupling 
between hydrodynamic modes in NE states considered 
here |55] . 

Since, in contrast to the Casimir pressure induced 
by NE temperature fluctuations, the Casimir pressure 
Pne induced by NE concentration exerts an actual force 
on the walls confining the liquid layer, we believe that 
the Casimir pressure induced by NE concentration fluc¬ 
tuations is a more promising candidate for an initial at¬ 
tempt to detect the phenomenon experimentally. As can 
be seen from Eq. Q, the effect can be enhanced by se¬ 
lection of a mixture with a small diffusion coefficient D 
and a large Soret coefficient St- This is the principal rea¬ 
son why in Tablej^pJ^g of l-methylnaphtalene-|-n-heptane 
with St = 1.73 x 10“^ [35] is much larger than 

in toluene-1-n-hexane with St = 0.32 x 10“^ K~^ [3(1] . 
The validity of linear non-equilibrium fluctuating dy¬ 
namics for the NE temperature and NE concentration 
fluctuations has been confirmed experimentally both by 
light scattering [Udiisa and by shadowgraph exper¬ 
iments m- Experimental evidence for the existence of 
a NE Casimir pressure would provide evidence for the 
validity of non-equilibrium nonlinear fluctuating hydro¬ 
dynamics. 

We note that similar NE concentration fluctuations 
and, hence NE Casimir forces, will also be present in 
liquid films with an isothermal concentration gradient or 
chemical-potential gradient [THl l42ti44] . Hence this kind 
of NE Casimir forces may be ubiquitous in nature. 
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